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Abstract 



The mixed cold-hot dark matter cosmological model (CHDM) with O,tot = 1 and a 
falling power-law initial spectrum of Gaussian adiabatic perturbations (n > 1) is tested 
using recent obserbational data. It is shown that its fit to the data becomes worse with 
the growth of n — 1, and may be considered as unreasonable for n > 1.1 for all possible 
values of the Hubble constant. Thus, the CHDM model with a falling initial spectrum 
is worse than the same model with the approximately flat (|n — 1| < 0.1) spectrum. 
On the other hand, the CHDM model provides a rather good fit to the data if n lies 
in the range (0.9 — 1.0), the Hubble constant Hq < 60 km/s/Mpc [Hq < 55 for n = 1) 
and the neutrino energy density Jlj^ < 0.25. So, the CHDM model provides the best 
possibility for the realization of the simplest variants of the inflationary scenario having 
the effective slope n ^ (0.95 — 0.97) between galaxy and horizon scales, including a 
modest contribution of primordial gravitational wave background to large-angle AT/T 
fluctuations of the cosmic microwave background (resulting in the increase of their total 
rms amplitude by (5 — 10)%) expected in some variants. A classification of cosmological 
models according to the number of fundamental parameters used to fit observational 
data is presented, too. 

1 Introduction 

It remains an ambitious goal of the inflationary scenario, as well as of any other fundamental 
cosmological theory of the early Universe, to explain all observed structure of the present- 
day Universe using a minimal number of additional microphysical "fundamental" constants, 
apart from those already khown from the particle physics. Of course, we don't know how 
many parameters is really needed to describe the whole Universe, so one can't say apriori 
that a theory having more parameters is worse than a theory with a less number of them. 
However, following the the Occam 's razor principle, classification of different cosmological 
models according to the number of additional phenomenological parameters used in them 
gives us a natural logical sequence of their consideration and comparison with observational 
data. We call these parameters fundamental, if they appear in basic equations (as in the 
inflationary scenario), not in initial conditions or other assumptions. 

As is well known, a power spectrum of perturbations producing the observed structure 
of the Universe is a product of an initial (primordial) spectrum formed in the early Universe 
and a transfer function C^(/c) which depends on the type of dark matter at present (e.g., on 
masses and concentrations of neutrinos). From the inflationary scenario point of view, the 
initial spectrum is completely determined by a phenomenological Lagrangian of an effective 
scalar fleld (or flelds) - inflaton(s) - at the de Sitter (inflationary) stage in the very early 
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Universe. Then parameters determining both the initial spectrum and the present dark 
matter content are fundamental and should be considered and counted on equal footing. If 
such a classification is applied to inflationary models (see, e.g., Starobinsky 1993), then a 
model of the first level having only one fundamental parameter - an amplitude of pertur- 
bations - is the CDM model with the approximately flat (Harrison-Zeldovich, or n ^ 1) 
spectrum of initial adiabatic perturbations. Because of theoretical considerations and ob- 
servational uncertainties, it is better to include "weakly-tilted" models with |n — 1| < 0.1 
into this class, too. Hereafter this model will be referred as the Standard Cold Dark Model 
(SCDM). There exists another model belonging to this level: the CDM model with the ap- 
proximately fiat spectrum of isocurvature fluctuations {n ^ —3). But that model has been 
known to be excluded by observations long ago, because it produces excessive large-scale 
AT/T fluctuations. Strictly speaking, SCDM has one more parameter which defines an 
amplitude of the approximately fiat spectrum of primordial gravitational waves and which 
is directly connected to the Hubble parameter H = a /a at the de Sitter stage (Starobinsky 
1979), a{t) being the scale factor of the Friedmann- Robertson- Walker cosmological model. 
However, this is a rather small effect which may be seen in a slight increase of large-scale 
AT/T fluctuations only (apart, of course, from a remote possibility of direct detection of this 
relic gravitational background), see the discussion section below. So, it is better to consider 
parameters connected with the gravitational wave background separately. 

At present, it is clear already that SCDM predictions, though being not far from obser- 
vational data (that is remarkable for a such a simple model with only one free parameter), 
still definitely do not agree with all of them. Namely, if the free parameter is chosen to fit 
the data on scales (100 — 1000)/i5o^ Mpc, discrepancy of about twice in perturbation am- 
plitude arises on scales (1 — lO)h'^Q Mpc, and vice versa {h^o — Hq/5Q, where Hq is the 
Hubble constant in km/s/Mpc). Thus, models of the next (second) level having one more 
additional constant have to be considered. Among these models, the best is certainly the 
mixed cold-|-hot dark matter model (CHDM), Shafi & Stecker 1984, for recent analysis see 
Pogosyan & Starobinsky 1993 (hereafter PS), Liddle & Lyth 1993, Klypin et al. 1993 and 
references therein. In this model, the hot component is assumed to be the most massive of 
3 neutrino species (presumably, r-neutrino) with the standard concentration following from 
the textbook Big Bang theory. Then the only new fundamental parameter is the neutrino 
mass m,y (masses of the other two types of neutrinos are supposed to be much less and, 
therefore, unimportant for cosmology). If, on the contrary, masses of two neutrino types 
are assumed to be comparable or even equal, the resulting model will belong to the third 
level, until the mass ratio will be either confirmed in laboratory experiments, or theoretically 
derived from some underlying theory (we shall return to the discussion of this case at the 
end of the paper), rriu is related to the energy density of the hot component (in terms of the 
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critical one) fli, by 

rriu = 23.31]^/i^o (1) 

for Try = 2.735K. The CDM model with the cosmological constant seems to be on the 
second place by a number of difficulties (and still marginally admissible) , and the two tilted 
CDM models with a power-law initial spectrum of adiabatic perturbations (with n < 1) 
and isocurvature ones (with n > — 3 to avoid excessive large-scale AT/T fluctuations) are 
marginally or completely excluded. 

Still the CHDM model with n ~ 1 is not without difficulties. The main of them is 
connected with later galaxy and quasar formation in this model as compared to the SCDM 
model. As a result, only a small region in the Hq — plain remains permitted (PS, see 
also a more pessimistic view in Cen & Ostriker 1994). Recently, this difficulty exacerbated 
due to the problem of producing sufficient number of damped Ly — a systems (Subramanian 
& Padmanabhan 1994, Mo & Miralda-Escude 1994, Kauffmann & Chariot 1994). However, 
latest analysis based on A'"-body simulations suggests that the latter problem may be solved 
if Qj, is taken smaller than it was supposed before: Q^, < 0.25 (Klypin et al. 1994) or even 
Oj, < 0.2 (Ma & Bertschinger 1994) for h^o = 1, in complete agreement with the restrictions 
on the model following from quasar and galaxy formation which were earlier obtained in PS 
using linear theory. 

Therefore, it is desirable to have more power on small scales in the CHDM model. This 
was one motivation for us to consider a CHDM model belonging to the next (third) level, 
i.e., having one more parameter. We assume the standard neutrino concentration and use 
one of three adjustable parameters as the neutrino mass as earlier. Then we are left with two 
parameters to characterize an initial perturbation spectrum. The most natural possibility 
is to assume a power-law spectrum of adiabatic perturbations with n ^ 1, then one of the 
parameters gives a rms amplitude of perturbations at some scale (say, at the present horizon 
scale), while the other defines the slope. For the reason stated above, we consider the case 
n > 1 in this paper (with some results relevant to the approximately fiat case |n — 1| < 0.1, 
too). The case n < 1 for the CHDM model has been already considered in detail in Liddle 
& Lyth 1993 and briefly mentioned in PS, with the conclusion that the "really" tilted case 
with n < 0.9 is excluded, but typical chaotic inflationary spectra with n ^ 0.95 — 0.97 (which 
we count as approximately flat ones) are possible. 

The other motivation to consider such a model is that the best ffi to the COBE data is 
given by n slightly larger than 1: n ^ 1.2 (Bennett et al. 1994, Gorski et al. 1994), although 
n = 1 lies inside la error bars (and, as advocated by Gorski et al. 1994, becomes even the 
best fit if the quadrupole is completely excluded), n > 1 is also needed to account for the 
results of the Tenerife experiment (Hancock et al. 1994). Of course, n significantly larger 
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than 1 produces too much power at small scales and, thus, may be rejected for different 
reasons. E.g., if the perturbation spectrum is assumed to maintain its power- law form up to 
very small scales where metric perturbations (gravitational potential) become comparable to 
unity, then the upper limit n < 1.4 follows from the consideration of production of primordial 
black holes (Carr, Gilbert & Lidsey 1994). Even if a falling power-law spectrum is assumed 
over a much smaller scale range of a few orders of magnitude, another upper limit n < 1.54 
follows from the absence of spectral distorsions of the cosmic microwave background (Hu, 
Scott & Silk 1994). For these reasons, we investigate the range n < 1.3 only. 

On the other hand, in spite of the advent of the Hubble telescope, there is still no general 
agreement about the value of the Hubble constant Hq. One methods produce the value of 
Ho around 50 km/s/Mpc (Branch & Miller 1993, Sandage et al. 1994), others lead to a 
significantly larger value (70 — 80) km/s/Mpc (Schmidt et al. 1994), while the accuracy 
of new methods based, e.g., on the Sunyaev-Zeldovich effect in rich clusters of galaxies, is 
still not enough to discriminate between these two cases (Birkinshaw & Hughes 1994). That 
is why we have to investigate a dependence of the model predictions on Hq additionally. 
We assume Hq in the range (40 — 80) km/s/Mpc. Of course, Hq is neither a new, nor a 
fundamental parameter, so it should not be counted in our classification scheme. 

Finally, to justify the use of the notion of fundamental parameters, we have to present 
at least one inflationary model producing a falling power-law spectrum over some range of 
scales. The simplest way is to take the well-known case of a test massive scalar field in 
the de Sitter background, that is equivalent to assuming the inflaton potential to be of the 
form y(0) = Vo + with ~ Hf = -C Mj, and the field being in the regime 

101 -C Mp [Mp = G~^/'^ is the Planck mass, and Ti — c — 1 is used in this paragraph). 
The inflationary stage in this model ends when drops to a very small value 0/ ^ Mp 
either as in the new inflationary model, or due to a second-order phase transition destroying 



Vq (Linde 1994). The exact expression for the slope is n — 1 = 3 — ^9 — 4^. The slope 

n = 1.1 corresponds to m = 0.38ifi (of course, then it is possible to use an approximate 
form of this formula: n ^ 1 + |^). Therefore, we may take m,,, ^ and as the three 
new fundamental parameters of the cosmological scenario considered. 
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2 CHDM Linear Perturbation Spectrum 



Following our previous notations for the case of the fiat initial spectrum, we write the linear 
power spectrum of perturbations of the gravitational potential $ in the form 

" 2^y^ (^) ■ ^CHBuik, t) , (2) 

specifying an initial spectrum produced at an early (inflationary) stage of the evolution of 

the Universe to be the power-law one: P<j,^in = 200^2^3 ("^) • The transfer function Cchdm 
was determined in our previous paper (PS) numerically by solving a system of the Einstein- 
Vlasov equations for the evolution of adiabatic perturbations with a neutrino component 
treated kinetically and a CDM component as dust. 

The resulting Cchdm depends on present values of the neutrino fractional density Q^, 
the Hubble constant Hq = bOh^o km/s/Mpc and the CMB temperature = 2.735t^ K 
as parameters, as well as on time t (redshift z). If we define Cchdm ^u,z) = Ccdm(^) • 
D{k, fli,, z), our numerical calculations can be described by the following fitting formula 

^^"^ - I 1 + {Bky - {Bkf + {BkY I ' ^-^^ 

with (5 — ^~^^4~^^" - ■ This formula satisfies asymptotic regimes for the transfer func- 
tion discussed in PS. Here ao/aeq — (1.681 Q^)"^ is the expansion factor of the Uni- 
verse from the matter-radiation equality moment teq until the present time to- The scale 
B = Rnr ■ (1 + Qt,o)/{fli, + fl^o), where Rnr = 10.80 h'r^ot'^ Mpc, is the neutrino non- 
relativization scale for Q^, = 1. For the critical value ^2,^ = ^2,^0 = 0.1435, neutrinos 
become nonrelativistic at t = tgq. The best fit of the only scale remained gives A — 

(1 + 10.912^^) - 0M65n^) I (1 + (9.2590^)^) with R^ = 69.52 hl^tl^ Mpc. The 

deviation of this fit from the calculated D{k) is better than 2% in the region k < 5/150 Mpc~^ 
for the present moment z ~ and fli, < 0.7, while staying within 5% for z < 30. We used 
the Bardeen et al. (1986) expression for Ccdm(^) in actual fitting. 

If the present CMB temperature is fixed (we use = 2.735 K), the model involved has 
four parameters, namely. A, n, fli, and Hq. If the adiabatic mode is solely responsible for 
the observed large-angle CMB fiuctuations, we have 

/Ary ^ A' r(3-n) r(/ + ^) 

\TJiJ 100 23-«r2(2-f)r(/ + 2-^^) ^' 
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(Bond & Efstathiou 1987). Note also a simple form of the corresponding large-angle corre- 
lation function (with the dipole components excluded) which we did not notice in literature: 



AT , . 



AT 



r(n)cos^^ 
2007r2 2'^-i(n- l)(2-n 

2 {?> — n){b — n){y — n 



a. 



10(n2 




l-n 



1— n 



2 (. ^ 3(n-l) ; 
1 + — — cos v 



3 — n 



3 — n 



(5-n) 

3(n- 1) 
(5 -n) 



cos-i? 



(5) 



-1 < n < 1. The 



where = -^{(^^^^ )• This expression is valid for both n > 1 and 
known correlation function for the n = 1 case (Starobinsky 1983) follows from here by 
limiting transition. If the tensor mode (gravitational waves) is responsible for a part of 

observed AT/T fluctuations, the value A derived from (^) for given ((^)^ ) serves as an 
upper limit on the amplitude of the adiabatic mode. However, we don't expect a noticeable 
tensor contribution for n > 1, unless a new parameter is introduced into the model that 
would shift it into the next, fourth level. 



3 Confrontation with Observational Tests 

Our way of comparison of the model with observational tests closely follows our previous 
paper (PS). However, here we confine our consideration to the following tests shown to be 
the most restrictive for the CHDM model in PS. 

1. The COBE measurement of large-angle AT/T fluctuations. 

2. Value of the total rms mass fluctuation as at R = 16/159^ Mpc (we use the index "8", 
not "16" , to be in accordance with standard notation caused by the habit of measuring 
Hq in units of 100 km/s/Mpc). 

3. The Stromlo-APM counts in cells (Loveday et al. 1992). 

4. Density of quasars (Efstathiou & Rees 1988, Haehnelt 1993). 

5. Large-scale peculiar velocities following from the POTENT reconstruction (Bertschinger 
et al. 1990, Dekel 1993). 
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Other tests included in the extended hst of PS do not lead to additional limitations on the 
model. 

We adopted the following numerical values for the observational data considered. 

The COBE result for the total rms value of the / = 4 multipole {AT/T)^ = (12.8 ± 2.3) 
liK/T^ (Wright et al. 1993) is used to put limits on the amplitude A, because it seems 
to be the most spectrum independent. For n = 1, this corresponds to the amplitude A = 
(4.38 ± 0.79) X 10~^ and the total quadrupole value Qrms-ps = (17.4 ± 3.1) /iK which are 
somewhat higher than those used in PS. 

The total rms mass fluctuation as = (^^(16/i^o^Mpc)^^ ^ is calculated for the lower limit 

of the COBE amplitude {AT/T)^ = 10.5 fxK/T^ (note that the inverse quantity cr^^ is usually 
understood as the biasing parameter b for optical galaxies). Based on cluster abundance data 
(White, Efstathiou & Frenk 1993), we consider ag < 0.67 as a conservative upper limit on 

The Stromplo-APM counts-in-cells represent nine data points for the mass variance in 
cubic cells cr^(/) over the range / = (20 — 150) h^Q Mpc. Fixing normalization by the con- 
dition (7^{l) = 1 for / = 25h^Q Mpc which agrees with the data, we eliminate a constant 
redshift correction (Kaiser 1987) and become able to compare directly the shape of the power 
spectrum with theoretical predictions for these scales. To find the best fit to the data, we 
formally applied the test with = 9 — 2 = 7 degrees of freedom for a fixed n. 

The standard model of quasar formation assumes that they result of formation 

of massive black holes in nuclei of galaxies with total masses (10^^ — 10^^) M0. Recent 
estimates of the fraction of mass in bound objects which can serve as quasar hosts at z = 4 
give /(> 10^^ M0) > 10~^ (Haehnelt 1993). In this paper, Haehnelt also presented the 
estimate for a fraction of mass in large galaxies: /(> lO^^Mg) > 10"^ at the same redshift 
z = 4. We use the simple Press-Schechter formalism to connect the mass fraction /(> M) 
to a linear mass fluctuation a{M, z) on a scale M at a redshift z: 

Then we get the limitation 

[(1 + z) a{M,z)]^^^>a6, (7) 

where a = 1.285 (1.132) for M = 10^^ (lO^^) ^g. The left-hand side of Eq. (g) depends 
on z only due to the ^-dependence of the transfer function There is no consensus on 
the threshold value 6c to be used. The standard one for top-hat fluctuations is 6c = 1.69, 
while Klypin et al. (1994) advocate for 6c = 1.4 as the best fit to a mass distribution in 
the CHDM model at high z. In the latter case, however, the Gaussian filtering was used to 
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calculate a{M) for a given mass scale. Let us note that the choice M = 10^^ M0 provides a 
tighter limitation than M = 10^^ M0. Therefore, we adopt with 6c = 1.4 (but with the 
Gaussian filtering) as a conservative restriction, while the limit on 10^^ Mq objects with the 
more standard 6c = 1.69 for top-hat fluctuations as, may be, more realistic one. Also, we 
found useful to use an amplitude-independent combination of the cxg limitation erg < a^: and 
the number of quasars test in the form 

[(1 + z) a{M, z)]^^, a-' > a6ca;' . (8) 

From the POTENT data, we use two values of large-scale bulk velocities v{80h'^Q Mpc) = 
(405 ± 60) km/s and v(120/i5o^ Mpc) = (340 ± 50) km/s (Dekel 1993). These values are in 
agreement with results obtained by other groups, see, e.g., Courteau et al. 1993 where the 
values w(80/i5o^ Mpc) = (385 ± 38) km/s and viUOh^f^ Mpc) = (360 ± 40) km/s for our bulk 
motion are presented. Note, however, that all these values refer to bulk flow velocities in 
our vicinity and, thus, may differ from their rms values calculated for the whole Universe 
("cosmic variance"). 

In Fig. 1, we display restrictions in the Hq — plane for several values of the slope of 
initial spectrum n = 0.85, 0.95, 1.1, 1.2 based on a combination of tests, namely the counts-in- 
cells cr^(/) values and the combined quasar density - ag condition (^. The case n = 1.0 was 
extensively discussed in PS (see Fig. 6 therein). The main conclusion made is that CHDM 
model parameters are restricted to the narrow range of a low Hubble constant Hq < 60 
km/s/Mpc and the neutrino fraction fij, = 0.17 — 0.28 for Hq = 50. The first of these limits 
reflects a problem with unavoidable high mass fluctuations at the IGh^o Mpc scale, as well 
as a wrong shape of the perturbation spectrum over the / = (20 - l50)h^Q Mpc interval if 
the Hubble constant is high. The upper bound on comes from the combined condition 
(H) that implies that the slope of the spectrum in the scale range {0.7h^Q^^ — IGh^^) Mpc 
cannot be too steep. The lower bound on Q^, arises both by matching of the spectrum slope 
to the counts-in-cells cr^(/) values (that shows more relative power on large scales than in 
the SCDM model), as well as from the result that cig is too high for low-Q,^ models with the 
COBE normalization. 

Fig. 1 shows how these results are affected by allowing an initial power-law spectrum to 
be non-flat: n 7^ 1. Two main conclusions can be drawn. The amplitude independent tests 
remain in essentialy the same mutual relation favouring somewhat higher values of Qi, for 
larger n as expected . For illustration, one may follow the point of intersection of the = 2 
contour with the limiting line for 10^^ Mq objects. This point moves from Hq = 62, Q^, = 0.3 
for n = 0.85 to Hq = 52, = 0.38 for n = 1.2. On the other hand, restrictions from the 
(Tg test become much tighter as the slope of initial spectrum increases. For n = 1.2, all 
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models with Hq > 50 are antibiased: 6 < 1. Thus, models with large n have too large mass 
fluctuations at the IGh^Q Mpc scale that cannot be compensated by increasing the neutrino 
fraction Q^y. 

We present another look at these results by plotting the n — Q^, cross-section of the 
parameter space for the same three tests in Fig. 2 . One can see that, for low values of 
the Hubble constant Hq = 40, 50, there is a significant degree of degeneracy between n and 
fl^, as far as the tests on the spectrum shape are considered, which do not disallow any 
value of n. Note, however, that if we adopt the more stringent quasar abundance condition 
/(> 10^^ M0) > 10^^ at 2; = 4, then we get an upper limit on Q^, for Hq = 50 roughly 
coinciding with the lower boundary of the best = 2 region for the counts-in-cells fit. So 
both tests are only in a marginal agreement for parameters close to the line SQ^, = n — 1/4. 

As we increase Hq, the fit to the counts-in-cells values begins to fail for large spectrum 
indexes, so that the best = 2 contour sets the limit n < 1.1 for Hq = 60, while for Hq = 70, 
this limit follows even from the more conservative condition ^ 7. Moreover, the best-fit 
contours tend to select higher values of the neutrino fraction than is allowed by the quasar 
test. In this way, for Hq > 55, no model can satisfy the condition (0) for M = 10^^ Mq and 
have < 2 (for Hq > 70, even < 7) simultaneously. If the less stringent quasar test 
for M = 10" M0 is used, it becomes possible to achieve the best fit for the counts-in-cells 
data for i^o = 60 if n < 0.9, but not for the Hubble constant as high as Hq = 70. On this 
basis, we conclude once again that the CHDM model is incompatible with high values of the 
Hubble constant Hq > 60 even for n 7^ 1 initial spectra. 

Although tests on the spectrum shape do set a tight limitation on the initial slope for 
high Hq, too, the strongest argument against high n comes directly from the erg condition. 
To achieve as < 0.67 for Hq = 50, one must restrict the model to n < 1, while for Hq = 60, 
even erg < 1 leads to n < 1.05. 

All the tests on the CHDM model become more in agreement with each other for Hq = AO. 
However, we are not sure that such low values of the Hubble constant are possible. Therefore, 
on the basis of Figs. 1, 2, we consider a small region of the parameter space in the vicinity 
of -f^o = 50, Qi, = 0.23, n = 0.95 as the best parameter set for the model. It corresponds to 
the neutrino mass rrii, ^ 5 eV. 

The Stromplo-APM counts-in-cells o"^(Z) are given in redshift space. Although we ex- 
cluded any constant redshift correction using these data as a test on the slope of the initial 
spectrum, it may be asked how the results would change if the redshift correction depended 
on scale significantly. Here we note that the counts-in-cells test in the form we used it serves 
primarily to indicate the necessity for relatively high Q^. Probably, the redshift correction 
can only increase with scale. Then underlying real space mass fluctuations depend less 
steeply on scale and are better fitted by lower Qi, models than the straightforward use of 
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cr^(/) predicts. 

To support the conclusions derived from the counts-in-cells test, we present contours 
of the direct fit of theoretical power spectra to a power spectrum reconstructed from the 
galaxy angular correlation function w{0) (Baugh & Efstathiou 1993, 1994) in Fig. 3. These 
data give the power spectrum directly in real space. We used their values for P{k) in the 
range k = (0.07 — 0.025)/i5o Mpc~^ where errors are not as large as on larger scales and 
nonlinear corrections are not yet as important as on smaller scales. In fact, we have only 4 
data points in the considered range which are fitted with the 3-parametric model (for fixed 
n or Hq). In Fig. 3, we compare contours of this one and the counts-in-cells tests in the 
Ho — Qi, plane for n = 0.95 and in the n — Qu plane for Hq = 50. First, we should note that 
no CHDM model fits the Bough & Efstathiou data too well (except for Hq = 40 or low n ). 
In particular, < 2 is achieved only for Hq < 50 if the spectral index is n = 0.95. Second, 

contours for these two tests have a rather close structural resemblance. Since the number 
of P{k) points used is small and the errors given by Bough & Efstathiou (1993) may be 
questioned, we don't use this test to set specific restrictions on the model. However, using 
a rather relaxed limitation < 4 for the angular correlation function test, we can confirm 
if not strengthen the lower limit on previously obtained from the counts-in-cells data for 
parameters in the most interesting region around Hq = 50, Q^, = 0.25. Other conclusions 
made as the deterioration of the fit with increasing of n or Hq are also confirmed by this 
test. 

In the previous consideration, we have not discussed the size of the allowed region in the 
amplitude A dimension of the parameter space. In Fig. 4, we present the Q^ — A cross-section 
of the parameter space with fixed i^o = 50. The shaded strip selects the adiabatic amplitude 
A following from the COBE result {AT/T)^ = (12.8 ± 2.3) /iK/T^. Now it depends on n 
(according to Eq. (0)) but not, practically, on Q^. This is also the case for the estimation of 
bulk velocities v{R) in spheres of radii R = 80, 120h'^Q Mpc. The condition ag < 0.67 and 
the quasar number test (0) for M = 10^^ Mq leave the triangular area left and below the 
point of intersection of solid curves as the allowed range of parametrs. Exactly this point of 
intersection produces the upper dashed line in Figs. 1,2. 

In Fig. 4, we see new aspects of the failure of the CHDM model to be successful for 
large n. Not only the Ug condition can't be fulfilled, but also large-scale bulk velocities 
become simultaneously too high to be compatible with ag < 0.67 for n > 1.1 and too low in 
comparison with the COBE amplitude for n > 1.2. 

We understand that there are two effects which might make these limits less strong. 
First, some part of the observed AT/T can be due to gravitational waves. Then the 
adiabatic amplitude A derived from the COBE data will be lower than that in Fig. 4. 
Second, bulk velocities in our vicinity can differ from average ones in the Universe by 
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cosmic variance. Therefore, assuming both that the velocities v{R) given by POTENT 
are at least 40% larger than their rms values and that the primordial gravitational wave 
background is responsible for the increase of the observed large- angle AT/T by at least 
1/3, one can, in principle, reconcile the amplitude A from these tests with ug < 0.67 for 
n = 1.2. However, the former assumption is equivalent to introducing one more funda- 
mental dimcnsionlcss parameter that shifts the model to the fourth level according to our 
classification. Really, since adiabatic and tensor fiuctuations are statistically independent, 
we have {{AT/Tf)tot = {{AT/Tf)^ + {{AT/Tf)g^. Therefore, to increase the rms value 
of {AT/T)tot by 40% or more, {{AT/Tf)^„ should be no less than {{AT/Tf)^ - a kind 
of additional "fine tuning" in the n > 1 case (in contrast to the n < 1 case where such 
a condition arises naturally). According to the main idea of classification of cosmological 
models presented in the Introduction, this does not mean that such a model is bad (because 
dimensionless fundamental parameters may be fine-tuned to some number) , it simply shows 
that more new significant fundamental parameters are necessary for explanation of observa- 
tional data than it was assumed initially. Thus, this way is not admissible if we want to stay 
among models of the third level. 

This can be illustrated using the simple infiationary model presented in the end of Intro- 
duction. Under the condition |0| <C Mp assumed earlier, the gravitational wave contribution 
is small: assuming n ~ 1 and using well-known expressions for the slow-roll motion, it is 
straightforward to show that 



for 2 ^ / < 30 (the numerical coefficient in round brackets is 6% more for / = 2, but 6% less 
for / = 3 and 9% less for / = 5, 6). Here the value of is, as usually, taken at the moment of 
the first horizon crossing during the inflationary stage (a; 60 e-folds before its end). Now let 
us try to get a larger value of T/S by assuming that |0| ~ Mp at this moment. To achieve 
this, we have to assume some specific relation between (f)f and Mp. Earlier, a value of the 
dimensionless constant (pf/Mp was not important for comparison with observational data, 
it was enough to assume that it was sufficiently small. But now this value becomes the new 
significant fundamental parameter of the model. Of course, then the power spectrum cannot 
be considered as an exactly power-law one. Note that, incidentally, this specific model of 
the fourth level does not achieve the aim of having T > S for n > 1.1, too, if we assume the 
abovementioned upper limit on the spectral index at very small scales following from the 
absence of excessive PBH formation: n < 1.4. This requires m < O.SHi, so, for scales where 
the local index n > 1.1, the temperature fiuctuation amplitude increases by \/l + ^ < 1-2. 
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For values of the spectral index n < 1, the situation with bulk velocities becomes inverse. 
For n < 0.85, they are higher than those following from the normalization to the large-scale 
l^jT. Now the possible effect of gravitational waves only worsens the discrepancy (for 
n = 0.85, in the framework of inflationary models with the exponential inflaton potential, 
the adiabatic amplitude A is 40% smaller than the one given in Fig. 4 due to this effect). 
Therefore, following Liddle & Lyth (1993), we exclude models with n < 0.9. 

4 Conclusions and discussion 

We compared the CHDM model with a falling power-law primordial spectrum of adiabatic 
perturbations {n > 1) with observational data. This model has one more adjustable dimen- 
sionless parameter (n — 1, which can be expressed in terms of an additional parameter of an 
inflaton potential) than the CHDM model with the approximately flat {n ^ 1) spectrum. It 
might be thought naively that the model with more parameters fits data better. Remarkably, 
we found just the opposite: fit to the data becomes worse with the growth of n — 1, and may 
be considered as unreasonable already for n > 1.1. 

Combining this with the previously known result that the CHDM model is hardly com- 
patible with observations for n < 0.9 (Liddle & Lyth 1993, PS), we arrive to the conclusion 
that the CHDM model requires the approximately fiat (|n — 1| < 0.1) primordial spectrum 
among all possible power-law spectra of adiabatic perturbations. This shows the robustness 
of the CHDM model with the simplest infiationary initial conditions. Our conclusion agrees 
qualitatively with that in the recent paper by Lyth & Liddle (1994). Moreover, the best fit 
to the data is achieved for n slightly less than 1, and around values of the effective slope 
expected in the simplest infiationary models: either with a scalar field with a polynomial 
potential V = n ^ 0.97 or F = 0.95 (chaotic infiation, Linde 1983), or in 

the higher-derivative gravity R-\- model (Starobinsky 1980) where n ^ 0.97, too, or in 
the new infiationary model with the Coleman- Weinberg potential (Linde 1982, Albrecht & 
Steinhardt 1982) with n ^ 0.95. 

Vice versa, the CHDM model presents the best possibility for the realization of these 
infiationary models, because the other alternative, the CDM-I-A model, has more serious 
problems. The most promising purely CDM models seem now to be based on a non-scale- 
invariant, step-like initial spectrum of adiabatic perturbations produced in more complicated 
infiationary models (Gootlobcr, Miickct & Starobinsky 1993, Peter, Polarski & Starobinsky 
1994), they belong to the third level in our classification. Thus, unexpectedly, the fate of 
the simplest infiationary models appears to be closely tied to the fate of the CHDM model. 

In addition, there is a place in the CHDM model with 0.9 < n < 1 for a modest but 
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noticeable gravitational wave contribution to large-angle — fluctuations expected for chaotic 
inflationary models with polynomial inflaton potentials (but not for the new inflationary 
model and the R + model). It leads to the increase in the total rms amplitude of these 
fluctuations by, e.g., \Jl + ^ — 1 ^ 10% for the V = ^ inflationary model (Starobinsky 
1985). Note that all the simplest inflationary models listed above do not require additional 
fundamental parameters to specify their predictions for (weakly scale dependent) values of 
the spectral index n — 1 and the ratio T/S. So, they belong to the second level if combined 
with the CHDM model for dark matter. The increase in ^ may be even found desirable in 
view of recent papers by Gorski et al. (1994), Banday et al. (1994) where further rise of the 
COBE results up to Qrms-ps = (19 — 20) fiK is suggested. On the other hand, it is clear 
that a significanly larger gravitational wave contribution is incompatible with the CHDM 
model. 

Of course, the most crucial confirmation of the CHDM model would be a direct or 
indirect (through neutrino oscillations parameters) discovery that the r-neutrino mass is 
really around 5 eV as predicted by the best fit to the model. But of no less importance is 
the precise determination of the Hubble constant because the model can't work for Hq > 60 
km/s/Mpc, and it is better to have its value around 50 km/s/Mpc. The third critical test 
for the model is the abundance of galaxies and quasars at large redshifts. 

Finally, let us return to the case of two (or even three) comparable neutrino masses 
mentioned in the Introduction. If neutrino concentrations are the standard ones and Q^, 
denotes the total neutrino energy density in terms of the critical one, then what stands in 
the left-hand side of Eq. (|l]) is actually the sum of masses of all neutrino species. In this 
case, the relative transfer function D{k, fl^, z) (|^) will have the same step-like form discussed 
in PS, but the transition region from one plateau to another will be shifted to larger scales 
due to increase of Rnr- So, in the first approximation, we may expect that the best fit is 
given by the condition that the sum of all neutrino masses should be about 5 eV. However, 
due to the abovementioned shift of characteristic scales in the transfer function, more caferul 
analysis of this higher level model is needed that will be carried elsewhere. 
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Figure captions 

Fig.l Restrictions in the Hq — parameter plane for n — 0.85, 0.95 (upper row) and n — 
1.1, 1.2 (lower row) following from: a) Stromlo-APM counts-in-cells test. Solid lines 
correspond to = 2, 7 contours, b) The cig condition. The values erg < 1, 0.67 are 
achieved left to the dashed (correspondingly right and left) lines, c) The combination 
of the as < 0.67 condition with quasar and galaxy formation conditions (dotted lines). 
The allowed region hes below the upper dotted hue if the hmitation set for objects of 
the mass M — 10^^ Mq is used and below the lower one if M = 10^^ Mq. 

Fig.2 Same tests as in Fig. 1 are given in the n — Q parameter plane for Hq — 40, 50 (upper 
row) and i^o = 60, 70 (lower row). 

Fig.3 Comparison of the contours for fits to the cloud-in-cells values (solid contours 
— 2, 7) and the power spectrum P{k) from the galaxy angular correlation function 
(dashed contours = 1, 2, 4). The left panel shows the Hq — Q^, plane for n — 0.95, 
the right one shows the n — Q^, plane for i^o = 50. 

Fig.4 Amplitude of perturbations A in the CHDM model as follows from potent data 
for bulk peculiar velocity v{R) [regions between dotted and dashed lines correspond 
to ±1(T error bars for v{80h'^Q Mpc) and f (120/i5o^ Mpc) respectively] compared to the 
COBE limits [shaded region]. The region above the solid line (i) is allowed by the 
quasar number condition /(> 10^^ M©, z = 4) > 10^"^. The region below the solid line 
(ii) is allowed by the condition erg < 0.67. 
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